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The  absorptivity  and  reflectivity  are  determined  from  the  Fresnel  coefficients  of  reflection  and  transmission  by  the 
complex  ray-tracing  method  for  a  plane  wave  incident  upon  the  surface  Z(x,  y)  of  an  isotropic  conducting  medium 
located  in  vacuum.  A  coordinate-dependent  real  refractive  index  of  the  form  AZ,  a,  u)  «  [qHZ,  a,  u)  +  2 q(Z.  a,  ui) 
cos  8in  +  l]17*  is  used,  where  q  and  cos  0i.  are  defined  in  part  I  of  this  series  [J.  Opt.  Soc.  Am.  A  5,  200  (1988)]. 
Numerical  results  are  presented  for  the  scattered  far  fields,  the  surface  waves,  and  the  scattering  and  absorption 
efficiencies  for  ellipsoidal  surfaces.  We  demonstrate  the  simplicity  and  accuracy  of  the  complex  ray-tracing 
method  over  the  scalar-potentials  method  in  solving  for  an  absorbing  medium  located  in  parallel  fields. 


1.  INTRODUCTION 


In  part  I  of  this  series1  we  generalized  the  expressions  for  q 
and  p  obtained  by  Stratton2  for  a  plane  wave  incident  upon  a 
plane  surface  to  the  general  case  of  an  arbitrary  surface.  In 
the  present  paper  we  apply  the  theory  developed  in  part  I  for 
ellipsoidal  surfaces.  Numerical  results  are  given  for  a  wide 
range  of  radiation  wavelengths  and  object  sizes  to  present 
quantitatively  the  different  physical  properties  of  the  scat¬ 
tered  Fields  resulting  from  the  irradiation  of  an  arbitrary 
metal  surface. 

The  paper  is  divided  into  five  sections.  In  Section  2  the 
mathematical  formalism  of  part  I  is  applied  to  ellipsoidal 
surfaces  located  in  a  parallel  field.  In  Section  3  numerical 
results  are  presented  for  the  local  and  integrated  absorp¬ 
tions.  In  Section  4  numerical  results  are  given  for  the  sur¬ 


A (Z,  <r,  u)  -  |0.5[(A2  +  B2),/2  -  A]|l/2. 


The  dependence  of  i»,  q.  A,  B,  'fr,  and  A  on  the  surface 
coordinates  Z,  x,  end  y;  the  conductivity  cr;  and  the  radiation 
frequency  u  are  assumed. 

An  ellipsoidal  surface  can  be  represented  by  a  coordinate 
function  of  the  form 


Z  *  (R2  -  y2y2  -  x2)1/2, 
where  y  is  the  eccentricity  and  is  given  by 


face  currents  and  the  scattered  far  fields.  Our  conclusions 


are  stated  in  Section  5. 


iv-  A.>  /  A.  /  ’  *  '  ■  r  \t 

2.  MATHEMATICAL  FORMULATION 


••V 

i  ►** . 


0<y  <=  1 


for  a  sphere  with  radius  R  and  with  its  center 
at  the  origin  of  the  Cartesian  frame  of 
coordinates 

for  an  infinite  cylinder  of  radius  R  with  its 
central  axis  coincident  with  the  y  axis, 
for  an  ellipsoid  with  principal  minor  axis  of  2 R, 
in  the  zx  plane,  and  principal  major  axis  of  26, 
on  the  y  axis  (6  *  R/y). 


The  letter  I  is  used  to  denote  equations  from  part  I.  Given  a 
plane  wave  with  a  propagation  constant  k  *  2r/\  traveling  in 
the  positive  z  direction  (0,  -  0)  and  using  Eq.  (I.IOa),  we  get 


The  geometries  of  the  three  surfaces  are  shown  in  Fig.  I.  To 
simplify  notation  we  define  the  following  parameters: 


q(Z,  a, «)  -  -($*  -  fA)/A  -  (l  +  Z2  +  Z2V 


The  real  part  of  the  complex  refractive  index  is  then  ob¬ 
tained  from  Eq.  (1. 15)  by  using  Fi  m  F-i  m  0  and  F3  «  I, 


*  yly,  g  -  (y2  -  Dy2,  u  •  (R2  +  gy2)1'2. 
The  surface  derivatives  are  then  given  by 
Zx  -  —x/Z, 


y(Z.  <r,  w)  -  (q2  +  2qcos9in  +  \)m. 


-w!Z, 


where,  from  Eq.  (I.3b),  cos  8m  ■  (I  +  Z,2  +  Zyr)~V2,  and  £,  f, 
♦,  ana  A  are  given  by  Eq.  (L8)  as 


A(Z,  <r,  w)  *  1  —  h2  sin2  8J£2  -  fW  +  (3a) 

B(Z,  e, «)  -  2k 2  sin2  8m  ff/({2  +  f)\  (3b) 

*(Z,  <r,  -  |0  5[(A2  +  B2)1'2  +  (3c) 


and  the  local  angle  of  incidence  is  defined  by 
cos  8-  ■  1/C  ■  Z/u, 


where  C  ■  (1  +  Z,2  +  Zy2) xn  [Eq.  (I.3b)j.  Substituting  Z„ 
Zy,  and  cos  0in  from  Eqs.  (6)  and  (7)  into  Eqs.  (1)— (3)  gives 
the  coordinate-dependent  refractive  index  for  the  surface 
Z(x,y).  Equations  (3a)  and  (3b)  give 


0740-3232/88/071105- 13*02.00 


©  1988  Optical  Society  of  America 


52 S 


r*y'*ir*rjrr  v*vv  vy  vv 


1106  J.  Opt.  Soc.  Am.  A/VoL  5,  No.  7/July  1988 


INCIDENT  PLANE 
WAVE 


P  (Xp.  Yp.  2p) 


INCIOENT  PLANE  WAVE 


\  !  \  * 
'•  ,'X 


- t 


/•x 


ILLUMINATEO 

SIDE 


NS 


A(Z,  a,  w)  *  1  -  A2[(x2  +  u/^/u2]^2  -  f 2)/({2  +  f2)2  (8a) 


5(2,  cr.  «)  -  2A2[(x2  +  yVu2]#/^  +  f2)2.  (8b) 

'F,  A,  7,  and  »  are  thus  determined  in  terms  A(Z,  a,  w)  and 
B(Z,  it,  w).  In  Fig.  2,  numerical  values  for  v  are  plotted  for 
aluminum,  using  R  *  0.1  cm;  y  m  0, 0.5, 1;  and  A  ■  0. 1  nm  and 
A  *  1.315  mh».  Figure  2(a)  shows  the  strong  coordinate 
dependence  of  v  at  short  wavelengths  (e.g.,  at  A  ■  0.1  nm  for 
an  aluminum  sphere  with  a  radius  of  0.1  cm  the  maximum/ 
minimum  ratio  for  »  is  1.0296).  Figures  2(b)-2(d)  show  the 
relatively  weak  coordinate  dependence  of  r  at  long  wave¬ 
lengths  (e.g.,  at  A  *  1.315  nm  the  maximum/minimum  ratio 
for  v  is  1.0010  for  the  above-described  aluminum  sphere). 

The  Fresnel  coefficients  of  reflection  and  transmission, 
written  in  terms  of  e,  are  as  follows*; 

R„  ■[(»*-  sin2  0jj  w  -  cos  8j[co*  ^(r2  -  sin2 
-  sin2  Jj/[(»2  -  sin2  9a)u2  +  cos  8 J 
X  ( cos  9m(r*  -  sin2  9iI)Ln  +  sin2  9  J,  (9a) 
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Fig.  1.  Geometry  of  the  three  ellipsoidal  metal  objects  used  to 
produce  the  numerical  results  in  Figs.  2-11.  (a)  An  infinite  cylin¬ 
der,  (b)  a  sphere,  and  (c)  an  ellipsoid  (7  »  0.5). 


R,  =  -[(k2  -  sin2  9Jm  -  cos  8j/ 

K*1  ~  sin2  8in)I/2  +  cos  8in],  (9b) 

Tp  »  2v  cos  8in/[(»2  -  sin2  8in)l/2 

+  cos  8iBj [cos  8m(v2  -  sin2  6  J 1/2  +  sin2  9  J,  (9c) 

and 

T,  -  2p  cos  9  J[( v2  -  sin2  9io)in  +  cos  8J.  (9d) 

The  reflectivity  and  the  absorptivity  obtained  from  the 
above  coefficients  are  as  follows: 

Pp  *  IK*1  -  sin2  9  Jin  -  cos  9  J 

X  [cos  8,„(r2  -  sin2  8J 1/2  -  sin2  9 J/ 

[(*-’  -  sin2  8to/‘ 1  -r  cos  9  J 

X  [cos  8iB(r2  -  sin2  9a)in  +  sin2  8J|2,  (9e) 

P.  ”  IK*1  -  »in2  8iB)u'1  -  cos  8  J/ 

[(**  —  sin2  9[a)in  +  cos  8JI1,  (90 


1 
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ap  -  4»2  cos  9a(r2  -  sin2  9J1'2/ 

[(/-sin2  cosSj2 

X  [cos  B  J.*1  -  sin2  9Jin  +  sin2  9  J2,  (9g) 

and 

a,  -  4  cos  0in(/  -  sin2  9m)in/ 

[(/  —  sin2  9ln)iri  +  cos  #in]2.  (9h) 

Again,  the  dependence  of  the  R,  T,  p,  and  a  terms  on  Z,  a, 
and  <*>  is  assumed.  The  local  intensity  of  the  energy  flow 
normal  to  the  surface  in  the  s  wave  (TE)  at  an  arbitrary  point 
on  Z  is  given  by  Eq.  (1.21)  and  for  the  p  wave  (TM)  by  Eq. 
(L22).  To  express  the  local  intensities  in  terms  of  the  sur- 
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face  geometry  we  start  with  Eqs.  (1.17)  and  (1.18)  and  substi¬ 
tute  $i  m  <Pi  m  0  and  0  for  the  polarization  angle  between  the 
electric  vector  on  the  incident  fields  and  the  x  axis.  This 
gives 

d  «  —cos  0»  +  sin  0 J,  (10a) 

t  -  (Zji  -  Zj)/(Z*  +  Z2)m,  (10b) 

t>  -  {ZJ+  zyj  +  (Z,2  +  Zy2)R\/[C(Z?  +  zyn  (10c) 
a  •  i  -  -(Z,  sin  4>  +  Zy  cos  0MZ,2  +  Z2)m,  (lOd) 

a  •  A  ■  (Zx  cos  0  -  Zy  sin  <t>)/C,  (lOe) 

and 

a  •  p  -  (-Z,  cos  0  +  Zy  sin  0)/[C(Zx2  +  Z  2)1/2].  (100 


i  I 

*  * 


Fif.  Z  Surface  coordinate  variation  of  the  real  refractive  index  r  for  aluminum  surfaces  with  R  »  0.1  cm  and  *  ■  0.  X  ■  0.1  * «m  in  (a),  and  X  • 
1.315  »in>  in  (b)-(d).  The  different  object  shapes  used  are  a  sphere  in  (a)  and  (c),  an  infinite  cylinder  in  (b),  and  an  ellipsoid  (7  -  0.5)  in  (d). 
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Fig.  4.  Variation  of  the  absorption  efficiency  with  the  radiation 
wavelength  for  spheres  of  aluminum,  copper,  iron,  and  gold,  calcu¬ 
lated  by  Eq.  (14). 

On  substituting  the  different  vector  products  from  Eqs.  (10) 
into  Eqs.  (1.21)  and  (1.22),  we  get 

j;  -  0.5(E,1/u)(ZI  sin  <P  +  Zy  cos  <P)2/[C(Z2  +  Zy2)}  (11a) 

and 

Jj  -  0ME*/v)(Zz  cos  0  -  Zy  sin  4>)V[C[ZX2  +  Z2)].  (lib) 

The  local  currents  on  the  surface  Z  are  given  by  Eqs.  (1.31) 
and  (1.32).  The  various  field  vectors  that  appear  in  the 
latter  two  equations  are  obtained  from  Eqs.  (1.19)  and  (10)  as 


1  «  ~[El(Zl  sin  <p  +  Zyco$  <P )/ 

(Z2  +  Z2)\(lZy-]Zx), 

(12a) 

*  -  -{(EJy))(Zx  sin  0  +  Zy  cos  0)/ 

(z,*  +  z;*)i(£z,+;zp. 

(12b) 

«  -[E,(Zt  coup-  Zy  sin  0)/ 

iz2 + z2mizx+ jzy). 

(12c) 

and 

V  -  [(Ei/vHZ,  cos  <f-Zy  sin  0)/ 

[Z2  +  Zy2)](lZ2-JZx).  (12d) 

Equations  (8M12)  fully  determine  both  absorption  and 
scattering  parameters  at  the  surface  Z,  as  shown  in  Sections 
3  and  4. 
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3.  ABSORPTION 

The  total  amount  of  power  absorbed  at  the  surface  Z  is 
obtained  by  integrating  the  local  absorbed  intensities  over 
the  surface  area  exposed  to  radiation,. 

P.  -  1 1  (apjj  +  atJ,')cos  M*dy.  (13) 

The  absorption  efficiency  is  then  defined  as 

7,dxdyY  (14) 

where  /,•  is  the  incident  intensity. 


(a) 


(b) 

Fig.  5.  Variation  of  the  Fresnel  reflection  coefficients  of  10.6-«tm 
radiation  at  the  surface  of  an  aluminum  sphere  with  R  «  0.1  cm  and 
#  »  0:  (a)  R,  and  (b)  Rr 
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Substituting  the  surface  derivatives  Zx  and  Zy  from  Eqs. 
(6)  and  the  absorptivities  from  Eqs.  (9g)  and  (9h)  into  Eqs. 
(11),  we  get 


J,'  •  §MZE?h\u)(x  sin  <t>  +  w  cos  <t>f/(x7  +  w2)  (15a) 


J.1  m  0.5(ZE*/riu)(,x  cos $  -  w  sin $)V(xJ  +  w1);  (15b) 


then  ths  total  absorbed  power  is  given  by 


P.  -  0.5(£,Vn)  (Z/u)V + 

Jn. 


X  (a,(x  sin  <t>  +  w  cos  $)* 


+  ap(x  cos  4  -  w  sin  $)*]dxdy. 


The  area  of  the  surface  Z  exposed  to  radiation  is  defined  by 
— 1  <  Hi  •  n  <  0.  Hence  Q«  is  determined. 

Numerical  values  for  the  local  absorption  (the  integrand 
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of  Eq.  ( 16))  are  plotted  in  Fig.  3  for  aluminum  surfaces.  The 
intensities  are  normalized  by  dividing  by  EiE’12 ij,  assuming 
that  Ej  is  independent  of  the  coordinates  for  a  plane-wave 
incidence.  In  Figs.  3(a)-3(d)  the  electric  field  vector  of  the 
incident  radiation  is  in  the  x  direction  (tf  *  0),  which  en¬ 
hances  the  absorption  along  that  direction.  In  Fig.  3(e)  the 
electric  field  points  in  the  y  direction  (<t>  m  90  deg),  and,  as 
expected,  the  enhanced  absorption  is  along  the  y  direction. 
By  examining  the  results  from  an  aluminum  sphere  with  a 
radius  of  0.1  cm  exposed  to  radiation  of  two  different  wave¬ 
lengths  [X  ■  1.315  Mm  in  Fig.  3(a)  and  X  -  10.6  m“  in  Fig. 
3(d)),  we  note  that  the  local  absorption  of  the  longer  wave¬ 
length  not  only  is  reduced  but  is  more  dependent  on  polar¬ 
ization.  This  observation  is  based  on  the  greater  elongation 
of  the  absorption  contours  of  the  10.6-Min  radiation. 

The  absorption  efficiency  )Eq.  (14))  is  computed  by  using 
an  automatically  adaptive  numerical  integration  method. 
The  results  are  plotted  in  Fig.  4  for  aluminum,  copper,  iron, 
and  gold  spheres  of  radius  R  m  0.1  cm.  The  absorption  cross 
section  Q.  is  obtained  by  multiplying  rja  by  the  geometrical 
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cross  section  of  the  object.  For  a  sphere  of  radius  0.1  cm,  Q„ 
ranges  from  0.03  cm2  at  A  «  0.1  nm  to  1.25  X  10-4  cm2  at  X  » 
0.1  cm.  Results  obtained  by  using  the  full-wave  treatment 
and  spherical  harmonics  (Ref.  3,  Fig.  13.15)  have  an  asymp¬ 
totic  value  for  Q„  of  0.0314  cm2  for  an  iron  sphere  of  radius 
0.1  cm.  This  result  agrees  reasonably  well  with  the  Qa  ob¬ 
tained  with  our  calculations.  The  damped  oscillation  (due 
to  diffraction)  of  the  curves  of  the  scattering  cross  section 
versus  the  radiation  wavelength  is  not  accounted  for  in  our 
method  of  solution,  which  is  based  on  a  complex  ray-tracing 
approach. 

4.  SCATTERING 

At  a  point  P(xp,  y^,  zp)  external  to  the  surface  2(x,  y),  the 
scattered  field  is  given  by  Eq.  (1.30).  Substituting  the  e*  and 
h‘  terms  from  Eqs.  (12)  into  Eqs.  (1.31)  and  (1.32)  and  using 


h„  *  [ifS  +  yj  +  Zpk)!D,  where  D  =  (xp2  +  yp2  +  zp2)1/2,  we 
get  the  following  expressions  for  the  surface  currents: 

A  X  e,  -  (Ei/C)(Sl{+  Sj  +  S3£)  .  (17a) 

and 

A  X  h,  -  (EjnQiTj  +  7V*  +  T3£).  (17b) 

The  S  and  the  T  terms  are  dimensionless  parameters  that 
account  for  the  physical  and  geometrical  properties  of  the 
surface,  the  degree  of  polarization,  and  the  radiation  wave¬ 
length.  The  S  and  T  terms  may  also  be  defined,  respective¬ 
ly,  as  the  partial  electric  and  magnetic  surface  waves  or 
currents.  They  are  obtained  from  Eqs.  (12),  (1.31),  and 
(1.32)  by 

=  [~ZX(  1  +  R,)(Zt  sin  0  +  Zy  cos  <t>) 

+  Zy(  1  -  Rp)(Zz  cos  0-2,  sin  0)]/(Zx2  +  Z,2),  (18a) 

Sj  -  l-Zy(  1  +  R,)(Zz  sin  0  +  Zy  cos  <t>) 

-  Zz{  1  -  RpHZ,  cos  ^  -  Zy  sin  0)]/(Zx2  +  Z,2),  (18b) 
S3  «  —(1  +  R,)(Zz  sin  0  +  Zy  cos  0),  (18c) 

^  -  [Zz(  1  +7y(Zx  cos  0  -  Zy  sin  0) 

+  2,(1  -  R,)(ZX  sin  0  +  Z,  cos  0)]/(Zx2  +  Z2),  (18d) 

Tj  -  [2,(1  +  Rp)(2x  cos  0-2,  sin  0) 

-ZJ1  -  R,)(Zz  sin  0  +  Z,  cos  <t>))l(Z2  +  Z2),  (18e) 

and 

7*3«  (l  +  Rp)(Zxcos0-Z,sin0).  (18f) 

Substituting  the  surface  derivatives  from  Eqs.  (6)  and  the 
surface  field  components  of  the  s  and  p  waves  from  Eqs.  (12) 
into  Eqs.  (18),  we  get  expressions  for  the  S  and  T  parameters 


in  terms  of  the  surface  coordinates: 

51  «  [—(1  +  R,){x2  sin  <&  +  xw  cos  4>) 

+  (1  —  Rp)(xw  cos  <t>  —  w2  sin  0)]/(x2  +  w2),  (19a) 

52  ■  [-(1  +  R,)(xw  sin  <t>  +  w2  cos  <t> ) 

-  (1  —  Rp){x2  cos  0  —  xw  sin  <t>)]/(x2  +  w2),  (19b) 

53  ■  (1  +  Rz)(x  sin  0  +  w  cos  0)/Z,  (19c) 

Tx  -  [(1  +  fip)(x*  cos  0  —  xw  sin  <t>) 

+  (1  —  R,)(xw  sin  <p  +  w2  cos  0)]/(x2  +  w2),  (19d) 

Tj  ■  [(1  +  Rp )(xu>  cos  $-w2  sin  <t>) 

-  (1  -  Jl,)(x2  sin  0  +  xw  cos  0)]/(x2  +  w2),  (19e) 

and 

Ts  ■  -(1  +  Rp)(x  cos  0  -  w  sin  0)/2.  (190 


It  then  follows  that  the  surface  currents  that  contribute  to 
the  scattered  fields  at  P  may  be  written  as 
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<b) 


«10' 


X  AXIS  *10" 

Fig.  7.  Contours  of  equal  scattered  intensity  at  the  surface  of  aluminum  objects.  $  •  0  and  X  «  1.315  pm  in  (a)-(c);  X  »  10.6  pm  in  (d).  The 
different  object  shapes  used  are  a  cylinder  in  (a);  a  sphere  in  (b)  and  (d);  and  an  ellipsoid  (y  -  0.5)  in  (c). 


X  (A  X  e.)  -  {Ej/DC)[t(y^S3  -  z^) 

+  y(ZpS,  -  x^Sj)  +  ^(XpSj  -  >^5,)],  (20a) 

Ap  X  [A,  X  (A  X  h,)l 

-  (Eji&c)  inVp^ 

-  (yp*  +  Zp*)T,  +  XpypTj 
+  J  [xpypTi  -  (Xpl  +  xp*)T,  +  ypXpTj] 

+  <  bp»pr2  -  (Xp*  +  V)r3  +  WJI-  <20b> 


The  three  components  of  the  scattered  electric  field  at  P(Xp, 
yp,  zp )  are  obtained  by  substituting  the  surface  currents  from 
Eqs.  (20)  into  Eq.  (1.30).  The  results  are 


EJP) 


-E\ik  exp(-iAD)/4*D3]  J  J  (Z/u)2|(ypSj 
-  ZpSj)D  -  [xpXpTj  -  (yp}  +  xp*)Tj  +  ipypTj| 


X  exp[tA(xxp  +yyp  +  Zzp)/D)  dxdy, 


(21a) 
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E„(P)  -  -E,[ik  exp(-ikD)/4*D3)  J  (Z/M)Jf(2pSl 

-  XpSJD  -  [xpypTl  -  (ip5  +  zp2)T2  +  ypJpTj | 

X  exp[ifc(xxp  +yyp  +  Z2p)/D]dxdy,  (21b) 

and 

£„(P)  »  -£;[ife  exp(-i/tD)/4irD3]  j  j  (Z/u)2|(xpS2 

-  ypSJD  -  [yp2pT2  -  (Xp2  +  yp2)T3  +  Vpr.)l 

X  exp(iA(xXp  +  yyp  +  Zzp)/D\dxdy.  (21c) 

Since  the  area  of  the  surface  Z  exposed  to  radiation  is  de¬ 
fined  by  —  1  <  hi  •  h  <  0  and  the  portion  of  this  area  that 
contributes  to  the  fields  at  P  is  defined  by  0  <  hp  ■  h  <  1,  fi, 
can  be  determined. 

The  depolarization  of  the  scattered  waves  is  determined 
by  defining  an  orthogonal  system  of  unit  vectors  along  the 
scattered  rays  similar  to  that  used  with  the  incident  rays 
[Eqs.  (10)].  Thus,  if  the  plane  of  incidence  of  the  scattered 
light  at  the  point  P  contains  the  unit  vectors  hp  and  £,  the 
orthogonal  unit  vectors  along  the  scattered  ray  at  P  are 

tp«itpX£/[l-</ip.£)2]l/2  (22a) 

and 

<2p  -  Ap  x  tp.  (22b) 

Substituting  fip  *  (xjl  +  yj  +  Zpi)/D  into  the  above  two 
equations,  we  get 

tp  -  (Zp/  -  y^)/D(zp2  +  y„2)in  (22c) 

and 

<?p  -  [-(Zp2  +  y2)l  +  Xpyp/  +  XpZp*]/D(Zp2  +  yp2)in.  (22d) 

Equations  (21)  and  (22)  give  the  angle  <t>p  that  the  electric 
field  vector  makes  with  the  plane  of  incidence  of  the  scat¬ 
tered  rays  at  P  as 

cos  0p  -  Re[-(Zp2  +  yp)E„  +  x^^,,  +  x pZp£„]/ 

i0(*p5+ypV/5iRe(£«)2 

+  RefEjy)2  +  Re<£J2H.  (23) 

The  scattered  wave  front  may  also  be  defined  by  surfaces  of 
constant  phase  shift  as  obtained  from  Eqs.  (21): 

+  -  tan-,[Ira(£„2  +  £./  +  £„V'5/ 

Re(£u*  +  £v2  +  £ul),/ll-  (24) 

The  above-described  treatment  may  also  be  extended  to 
the  case  of  an  aplanar  light  beam  by  taking  £,  as  a  function  in 
the  radial  coordinate.  For  example,  for  a  Gaussian  beam,  £, 
assumes  the  form  £g  exp[-(*2  +  y2)/wo2]. 


The  following  observations  of  the  local  and  far  scattering 
can  now  be  made: 

(a)  In  Fig.  5,  in  which  the  Fresnel  reflection  coefficients 
R,  and  Rp  [Eqs.  (9a)  and  (9b))  are  plotted,  we  can  see  that  the 
reflection  of  the  local  p  waves  (TM)  increases  at  the  central 
parts  of  the  surface  while  that  of  the  s  waves  (TE)  increases 
on  the  periphery. 

(b)  In  Fig.  6,  the  partial  surface  waves  [Eq.  (17)]  are  of 
the  primary  and  the  secondary  orders.  Partial  waves  of  the 
third  and  higher  orders,  obtained  by  Mie  (Ref.  3,  Fig.  13.9), 
have  no  correspondence  in  the  ray-tracing  approach. 

(c)  The  normalized  local  fraction  of  the  energy  scattered 
at  the  surface  is  obtained  by  an  equation  similar  to  Eq.  (13) 
but  with  reflectivities  p  in  the  place  of  the  absorptivities  a. 
The  numerical  results  are  plotted  in  Fig.  7  for  three  alumi¬ 
num  surfaces.  Comparing  the  dependencies  of  the  scatter¬ 
ing  contours,  such  as  that  shown  in  Fig.  7(a),  with  the  ab¬ 
sorption  contours  [Fig.  3(a)],  we  find  that  absorption  con¬ 
tours  are  more  dependent  on  the  polarization  and  more 
elongated  along  the  direction  of  polarization. 

(d)  The  integrated  scattering  efficiency,  defined  by  an 
expression  similar  to  Eq.  (14)  but  with  the  scattered  power 
in  the  place  of  the  absorbed  power,  is  plotted  in  Fig.  8.  At 
wavelengths  longer  than  4  pm  it  is,  to  a  good  approximation, 
practical  to  assume  that  the  four  metals  considered  are  infi¬ 
nitely  conducting  media  with  near-unity  scattering  efficien¬ 
cies. 

(e)  The  scattered  far-field  intensity  is  plotted  in  Fig.  9. 
Figures  9(a)-9(c)  show  the  scattered  power  from  an  alumi¬ 
num  ellipsoid  illuminated  by  iodine  laser  light  of  a  unit 
intensity  (A  *  1.315  Mm).  At  this  wavelength,  the  maximum 
scattered  power  (as  high  as  70%  of  the  incident  intensity)  is 


x ;  IN  pm 


Fig.  8.  Variation  of  the  scattering  efficiency  with  the  radiation 
wavelength  for  sphere*  of  aluminum,  copper,  iron,  and  gold. 
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Fig.  9.  Scattered  far-field  intensity  from  aluminum  surfaces  with  R  ”  0.1  cm  and  $  *  0:  X  »  1.315  >im  in  (a)-(c)  and  in  (g)-(i);  X  «  10.6  nm  in 
(d)-(f);  X  ■  1  cm  in  (j)-(l).  The  target  geometries  are  an  ellipsoid  (y  “  0.5)  in  (a)-(f)  and  a  sphere  m  (g)-(l). 


contained  in  two  sidelobes  in  the  zy  plane.  Maximum  back- 
scattering  occurs  at  25%  in  the  vicinity  of  the  surface.  Fig¬ 
ures  9(d)-9(f)  represent  results  for  COj  laser  radiation  (X  =* 
10.6  *rin).  Here,  the  amplitude  of  the  scattered  power  is 
much  less  (maximum  reflected  intensity  3  X  10-3)  than  for 
that  of  the  iodine  laser,  but  the  side  scattering  perpendicular 
to  the  plane  of  polarization  is  still  an  important  feature. 
Figures  9(g)— 9(1)  show  the  intensity  of  the  far-field  scattered 
radiation  from  aluminum  spheres  irradiated  by  the  iodine 
laser  [Figs.  9(g)— 9(i)]  and  by  the  microwave  radiation  [Figs. 
9(j)-9(l)[. 

(f)  Figure  10  shows  contours  of  equal  intensity  of  the 
microwave  radiation  scattered  from  aluminum  spheres  and 
is  intended  to  demonstrate  the  symmetrical  properties  of  the 
fields.  Contours  are  plotted  in  each  of  the  three  Cartesian 
planes  to  show  the  symmetry  of  the  scattered  fields  around 
H,  the  direction  of  propagation;  d,  the  direction  of  polariza¬ 
tion;  and  the  surface  orientation.  The  scattering  contours 
are  slightly  distorted  from  circular  curves  in  the  xy  plane 
and  elliptical  curves  in  the  yz  plane,  and  the  maximum 
scattering  intensity  is  as  low  as  10~8  of  the  incident  intensity. 

(g)  Figure  11  is  the  result  of  rotating  the  plots  in  Figs. 
9(i)  and  9(k),  respectively,  around  the  x  axis  and  the  y  axis  to 
show  the  forward  scattering  patterns.  The  dip  in  the  scat¬ 
tered  power  intensity  in  the  zx  plane  [Fig.  11(a)]  is  due  to  the 
enhanced  absorption  in  that  plane  (the  incident  electric 
vector  points  in  the  x  direction).  In  Ref.  4,  Fig.  8,  the  dip  in 


the  scattered  power  occurred  in  the  zy  plane  and  was  de¬ 
scribed  by  the  diffraction  effect,  since  the  sphere  used  there 
had  dimensions  that  were  comparable  with  the  radiation 
wavelength.  Also,  in  Ref.  4,  Fig.  4,  with  wo  -  4.18603X  and  R 
m  0.2093A,  the  beam  spot  size  ui0  was  taken  to  be  large 
enough  to  cover  the  whole  sphere,  and  the  maximum  scat¬ 
tered  power  occurred  in  the  zy  plane.  The  results  from  Ref. 
4  are  similar  to  those  obtained  for  a  plane  wave,  plotted  in 
Figs.  9(c)  and  9(f). 

5.  CONCLUSION 

Based  on  a  complex  ray-tracing  approach,  we  obtain  concise 
analytical  expressions  for  the  far  scattering  fields  that  ex¬ 
plicitly  include  the  surface  function  Z(x,  y)  and  its  physical 
properties,  Eqs.  (19)  and  (21).  We  have  demonstrated  the 
efficiency  of  this  approach  in  solving  for  the  absorption  and 
scattering  parameters  of  large  conducting  surfaces.  Accu¬ 
rate  quantitative  descriptions  of  the  magnitude,  phase,  and 
depolarization  of  the  scattered  far  fields  and  the  surface 
currents  are  presented.  The  results  obtained  are  beneficial 
in  the  design  of  metal  reflectors  and  in  the  study  of  the 
interaction  of  electromagnetic  fields  with  complex  conduct¬ 
ing  surfaces,  as  in  the  case  of  depolarization  of  the  backscat- 
tering  that  is  due  to  surface  roughness.  The  proposed  anal¬ 
ysis  can  be  .xtended  easily  to  the  case  of  aplanar  laser  beams 
such  as  Gaussian  and  higher-order  mode  beams. 
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Fig.  10.  Contours  of  equal  scattered  far-field  intensity  in  the  three 
Cartesian  planes  from  an  aluminum  sphere  with  /?•»  0.1  cm,  <t>  »  0, 
and  \  «  1  cm.  The  numbers  on  the  contours  are  the  logarithmic 
values  of  the  scattered  intensities  for  a  unit  intensity  of  incident 
light,  (a)  xz  plane;  (b)  yz  plane;  (c)  xy  plane. 
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